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ABSTRACT 

We consider the problem of self-regulated heating and cooling in galaxy clusters and the im- 
plications for cluster magnetic fields and turbulence. Viscous heating of a weakly collisional 
magnetised plasma is regulated by the pressure anisotropy with respect to the local direction of 
the magnetic field. The intracluster medium is a high-beta plasma, where pressure anisotropics 
caused by the turbulent stresses and the consequent local changes in the magnetic field will 
trigger very fast microscale instabilities. We argue that the net effect of these instabilities will 
be to pin the pressure anisotropics at a marginal level, controlled by the plasma beta parame- 
ter. This gives rise to local heating rates that turn out to be comparable to the radiative cooling 
rates. Furthermore, we show that a balance between this heating and Bremsstrahlung cooling 
is thermally stable, unlike the often conjectured balance between cooling and thermal conduc- 
tion. Given a sufficient (and probably self -regulating) supply of turbulent power, this provides 
a physical mechanism for mitigating cooling flows and preventing cluster core collapse. For 
observed density and temperature profiles, the assumed balance of viscous heating and radia- 
tive cooling allows us to predict magnetic-field strengths, turbulent velocities and turbulence 
scales as functions of distance from the centre. Specific predictions and comparisons with 
observations are given for several different clusters. Our predictions can be further tested by 
future observations of cluster magnetic fields and turbulent velocities. 

Key words: galaxies: clusters: intracluster medium - magnetic fields - instabilities - turbu- 
lence - diffusion - conduction 



1 INTRODUCTION 

Early X-ray observations of galaxy clusters and the intracluster 
medium (ICM) indicated radia tive losses large enough to lead to 
cooling flows (see reviews by Sarazin 1986, 1988). Mass deposition 
rates Mqf due to presumed cooling flows were estimated to be as 
much as ~ 10^ Mp; yr~^ in some clusters |Cowie & Binnevll977l : 
iFabian & Nulsenll 19771 ; iMathews & Bregmanlll978h . The cooling- 
flow model also predicted copious iron line emission from tempera- 
tures between 10^ and 10^ K. However, at the time there was little 
direct evidence for mas s dropout in a ny spectral band other than 
X-rays (for a review, see lFabianlll994) . 

More recent high spectral resolution X-ray observations failed 
to detect the expected iron-line emission and constrained the 
centr al temperature to be ~ 1/3 of the bul k cluster tempera- 
ture jfor a review, see Peterson & Fabiarj|2006h . The spectroscop- 
ically determined ma ss deposition rates are < 0.1 A/cf (e.g. 
IVoigt & FabTai] l2004h . This is despite the fact that the cooling 
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time at r < 100 kpc is less than a Hubble time in > 70% 
of clusters (e.g. E dge. Stewart & Fabian 1992; Per es et al. 1 1 19981 ; 
Sand erson. Ponman & O' SuUivanl |2006| ; IVikhlinin et al.l l2006h . 
This discrepancy is the so-called 'cooling flow problem.' 

Some heating mechanism must therefore be balancing radia- 
tive cooling in the 'cool-core' clusters. A wide variety of heating 
and heat transport schemes have been considered, including ther- 
mal energy from the outer regions of the clu ster be ing transpo rted to 

.iBinnev & Cowi3ll98ll; 



the central cooling gas by conduction (e.g. 



i Tucker & Rosner 1983'; Bertschinger & Meiksen 1986; Voigt et all 



(Fabian, Voigt & Morris 2002; Zakamska & Narayai 
energy in jets, bubbles, cosmic rays, outflows and/or 
radiation from a central active galactic nu cleus (AGN) either 
resulting in turbulent d iffusion of heat (e.g. lOeiss & JustI 1 199^ ; 
iKim & Naravai] l2003j) and/or being thermalised via dissipa- 
tion of tur bulent motions, sound waves and/or gravitational 
modes (e.g. Loewenstein & FabianI 1990; Binney & Tabor 199^;^ 



Ciotti & Ostrike^ '2001' 



Ruszkowski, Brti ggen & BegelmanI 
Voit & Donahue 2005; Fab ian et alj |2005|; 
Brtiggen. Ruszkowski & HallmanI 120051 ; IChandran & Rasera 



2004allbl; 
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I2OO7I ; iMcCarthv et al. 2008b: dynamical friction from galaxy 
wakes (e.g. ISchiDpej|l974 iLea & De Young|ll976l : lMillej|l986l : 
Just et al.' 199(]^ or som e combin ation of these ( e.g. Briiggen 
2003 : Dennis & Chandran 2005: C onroy & Ostriker 2008). Much 
of recent work has focused on either thermal conduction, AGN 
heating, or both. 

Thermal conduction alone cannot be the solution to the cool- 
ing flow problem across the full range of masses due to its 
steep temperature dependence jVoigt & Fabianl2004l : lKaastra et al.l 
[2004: Pope et al. 2006). Even in hot systems where conduction 
is potent, fine tuning of the suppression factor / (the fraction by 
which Spitzer conductivity i s reduced due to, e.g., ta ngled mag- 
netic field lines) is required teregman & Davidlll988[) . Moreover, 
if the thermal conduction has the same temperature dependence as 
the Spitzer conductivity (i.e. if / is a constant) for a given ICM 
atmosphere, the resulting equiUbria are thermally unstable (e.g. 
iBregman & D"avidlll988l : 1Sokeill2003l : Ixim & Naravanll2003bl : see 
further discussion in S ec tions [2 . 3 1 and IT6I 1 . 

More recently, however, there has been renewed interest in the 
possibility that thermal conduction may provide sufficient heating 
to stably counteract the effects of radiative cooling. This has gone 
hand-in-hand with a dramatic increase in our understanding of 'di- 
lute' (i.e. only weakly collisional) plasmas, due to the apprecia- 
tion that even very weak magnetic fields introduce an anisotropy 
into heat fluxes. One important consequence is that the criterion for 
convective instability changes t o one of temperatu r e, rather than en- 
tropy, increasing downwards ( B albu sl2OO0l I2OO ih . lOuataertI ( |2008|) 
generalised Balbus's (2000) analysis and found that a heat-flux 
buoyancy-driven instability (HBI) occurs for upwardly-increasing 
temperature profiles as well, so long as the magnetic field is not en- 
tirely horizontal (orthogonal to gravity). Balbus & Reynolds ( 200^ 
conjectured that the nonlinear HBI is self-regulating and drives a 
reverse convective thermal flux, both of which may mediate the sta- 
bilisation of cooling cores. Numerical simulations of the HBI have 
been performed by Parrish & Ouataert' ('2008'), Bogdan ovic et al.l 
12009) and Parrish, Ouataert & Sharma (2009) with applications to 
the ICM. It was discovered that the HBI acts to rapidly reorient field 
lines to insulate the core, undermining the role of thermal conduc- 
tion and causing a cooling catastrophe to occur. Subsequent work 
has shown that a moderate amount of turbulent driving may help 
regulate the HBI and allow thermal channels to remain open, po- 
tentially stabilising the core against collapse jSharma et al.l l2009l : 
IParrish Ouataert & Sharmall201(]| : lRuszkowski & OlJl2010l) . 

There are several reasons to believe that AGNs play an impor- 
tant role in regulating cooling. In many clusters, the AGN energy 
output inferred from radio-emitting plasma outflows and cavities 
is sim ilar to the X-ray cooling rate of the central gas ( Fabian et al.l 
2000: Kanov. Saraz in & Hicks|[2006l : iMcNamara & NulsenI 120071 : 
Forman et al. .2007.) . Moreover, > 70% of cool-core clusters 
harbour radio sources at their centres, w hile < 25% of clus- 
ters without cool cores are radio loud dBurnj Il990h . provid- 
ing strong circumstantial evidence for a connection between the 
processes that fuel the radio emission (such as AGNs) and the 
X-ray emission from the cooling gas. Models of self-regulated 
heating from AGN have been constructed (e.g. Ciotti & OstrikeJ 



2001: Ruszkowski & Begelman 2002: Brighenti & Mathews"2003T 
Kaiser & Binney 2003: Hoeft & Bruggen 2004: Guo & Oh 2008: 
Briiggen & Scannapiecoll2009l) in which AGN activity is triggered 
by cooling-induced gas accretion toward cluster centres, increas- 
ing AGN heating and halting the collapse. Episodic outflows from 
AGN are thought not only to quench cooling and condensation in 



clusters, but also to limit the maximum luminosit y of galaxies a nd 
regulate the growth of black holes at their centres ( lBinney||2005[) . 

While AGN activity is fundamentally linked with the observed 
prese n ce of radio bubbles an d/or X-ray cavities (e.g. Birzan et al.l 
I2OO4I: iDunn & Fabian! |2006|) . it is currently unclear how the 
AGN energy is actual ly thermalised (e.g. see the introduction of 
IVoit & DonS^l2005l for a review of possibilities). This question 
can only be answered once knowledge of the effective viscosity 
of the ICM is acquired. The ICM hosts subsonic turbulence and 
magnetic fields with energy density comparable to that of the mo- 
tions. Bot h of these should affect the vi scosity of the ICM (see 
review bv lSchekochihin & Cowievll2006h . In particular, the pres- 
ence of a magnetic field alters the form of the viscosity when the 
ratio of the ion cyclotron an d collision frequencies is much greater 
than unity teraginskiilll965h , a condition amply satisfied in galaxy 
clusters. As a result, the transport properties of the ICM become 
strongly dependent on both the geometry and strength of the mag- 
netic field, as well as on microscale plasma instabilities that are 
likely to occur ubiquitousl y in the ICM (e.g. firehose and mirror : 
Schekochihin et al. 2005: Lvutikovl2007l:ISchekochihin et al.l200^ 
l201Ql : iRosin et al.li2010l) . 

In this paper, we investigate the effect these plasma effects 
might have on the large-scale transport properties of the ICM. 
Specifically, we argue that parallel viscous heating, due to the 
anisotropic damping of turbulent motions, is regulated by the satu- 
ration of microscale plasma instabilities (e.g., firehose and mirror) 
and can balance radiative cooling in the cool cores of galaxy clus- 
ters in such a way as to ensure thermal stability (Section|2]l. Given 
observed densities and temperatures, this balance implies specific 
values for central magnetic-field strengths and radial profiles of the 
rms magnetic field that are in good agreement with current ob- 
servational estimates and that lend themselves to testing by future 
observations (Section|3]l. We also show that, under the reasonable 
assumption that turbulent kinetic and magnetic energies are com- 
parable to one another, cluster profiles for the turbulent velocity 
and the characteristic turbulence scale may be derived. The spe- 
cific case of A 1835 is considered as a typical example in Section 
13.41 Since the fundamental plasma-physical processes we appeal to 
are in principle universal in both cool-core and non-cool-core (i.e. 
unrelaxed) clusters, we also calculate in Section [33] the magnetic- 
field strengths and turbulence characteristics for some non-cool- 
core clusters. They turn out to be in good agreement with current 
observational estimates. Thermal conduction and the robustness of 
our results with respect to it are briefly discussed in Section [T6l 
In Section [4] we close with a discussion of our results and their 
limitations. 



2 COOLING AND HEATING OF THE ICM 
2.1 Radiative cooling 

We follow iTozzi & Norin^] ( l200lh in approximating the radia- 
tive cooling rate (per unit volume) of the ICM (as determined by 
[Sutherland &Do pita 1993.) by 



Q = nincA(r), 

where the cooling function is 

A(r) = 10"^^ erg s"^ cm^ 
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1 keV 
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Our notation is standard: ni (We) is tlie number density of ions 
(electrons) and T is tlie temperature (in energy units through- 
out the paper). We also define, for future use, the total number 
density n = ni + rig. The numerical constants Ci = 0.086, 
C2 ~ 0.58 and C3 = 0.63 are selected to correspond to an 
average metallicity Z — 0.3 Zq, for which the mean mass per 
particle is jj, — 0.597mp, the mean mass per electron is /ie = 
uinJne) = l.lSOmp and ni = 0.927nc (Sutherland & Dopital 
Il993h . The cooling function is dominated by Bremsstrahlung above 
r ~ 1 keV and by metal lines below T ~ 1 keV. 

Since the temperature equilibration time between ions and 
electrons fi_c,cq ~ 10 kyr near the centres and ~ 1 Myr near the 
temperature maximum of cool-core clusters, which is smaller than 
all other timescales that will be relevant to us (see Section [T6l >. we 
assume Ti — Tc — T. Even in unrelaxed clusters like Coma, where 
the te mperature is relatively high (T ~ 8.2 keV; lArnaud et al.l 
l200lh and the electron density is relatively low (rie ~ 3 x 10 



' 4 X 10"^ cm"^; see Section[331l, ti- 



2-170 Myr. 



Normalised to conditions representative of the centres of cool- 
core clusters, the radiative cooling rate (per unit volume) is 



Q~ 



1.4 X 10"^^ 



( V 

VO.l cm-3/ 



T 



2keV 



1/2 



erg s cm 



(3) 



in the Bremsstrahlung regime (T > 1 keV). 



2.2 Parallel viscous heating 

There is a rapidly growing body of observational evidence for the 
presence of appreciable magnetic fields in the ICM (for a review, 
see [Carilli & Tavlor 2002). Randomly tangled magnetic fields with 
strength _B ~ 1-10 /iG and characteristic scale ~ 1 - 10 kpc 
are consistently found, with fields in the cool cores of cooling-flow 
clusters somewhat stronger than elsewhere. 

The presence of a magnetic field alters the form of the ther- 
mal pressure when Q.\/vu 2> 1, where f^i = eB/m\c is the 
ion cyclotron frequency and = 'i^/Tini\is_e^ /?,m\^'^T^^'^ is the 
ion-i on collision frequency; An is the ion-ion Coulomb logarithm 
tera ginskii 1965). This is certainly the case in the cool cores of 
galaxy clusters, where typical values of the electron density Jio, 
temperature T and magnetic-field strength B imply 



Q.i 10 f B 



0.1 cm- 



T 



2keV 



3/2 



(4) 



As a result, thermal pressure becomes anisotropic with respect to 
the local magnetic field direction b: 



P = Px (I - bb) +P||bb = p\ + er, 



(5) 



where p± (p||) is the thermal pressure perpendicular (parallel) to 
the local magnetic field, p = {2/3)p± + (l/3)p|| is the total ther- 
mal pressure, I is the unit dyadic and we have defined the collisional 
viscous stress tensor 



a^-{bb--\ 



(6) 



In the Braginskii (i.e. collisional) limit, appropriate for the large- 
scale motions in the ICM since their dynamical timescales are 3> 
i'i7^ S> the ion contribution to the viscous stress dominates 
that of the electrons by a factor proportional to (mi/mc)'^'^^. Thus, 
in what follows, we neglect the electron contribution to the viscous 
stress. 



The viscous stress tensor appears both in the momentum equa- 



tion, 

du _ 
mmi— = - V- 



BB 

4-K 



(7) 



as a form of momentum transport, and in the energy equation, 

3 dr 

^n— = -nTV-n - miVu - V-q^ - nincA(r), (8) 

as a form of heating. In these two equations, g is the gravitational 
acceleration, is the electron collisional heat flux and d/dt = 
d/dt + u-'V is the convective derivative. We assume an ideal gas 
equation of state, p — nT (both for each species and for the two 
combined because Ti = Tc — T). The electron contribution to 
the collisiona l heat flux domin ates that of the ions by a factor of 
(mi/mc)^/^ ( lBraginskiilll965h . 

Differences between the perpendicular and parallel pressures 
in a magnetised plasma are due to the conservation of the first adi- 
abatic invariant for each particle, /i = mv\/2B — const (on 
time scales ^ i^f^). Therefore, any change in the field strength 
must be accompanied by a corresponding change in the perpen- 
dicular pressure, p±/B ~ const. In a turbulent plasma such as 
the ICM, time-dependent fluctuations in the magnetic-field strength 
are inevitable. Accordingly, a patchwork of regions of positive or 
negative pressure anisotropy will emerge, corresponding to locally 
increasing or decreasing magnetic-field strength. If the pressure 
anisotropy |px— P||| ^ B'^ /A-k, firehose and mir ror instabilities are 
triggered at spatial and temporal microscales ( Schekochihin et all 
2005, and references therein) Q Equations Q and ((8) break down 
at these scales, and the perpendicular and parallel pressures must be 
determined by a k inetic calculation (e.g. .Schek ochihin et a l. 201o|; 
iRosin et alfcoiOh . 

It is usually the case that the pressure anisotropy - and thus 
the viscous stress - is regulated by the nonlinear evolution of these 
microscale instabilities, which t end to pin t he pressure anisotropy 
at marginal stability values (see iRosin et alji2010i . and references 
therein): 



Ai = 



-P\\ 



Pi 



A' 



(9) 



where ^ = — 1 for the firehose instability or 1/2 for the mirror 
instability0 The (ion) plasma beta parameter 



ft 



STTOiT 



75 



B 



10 AtG 



0.1 cm- 



T 



2keV 



(10) 



is the ratio of the (ion) thermal and magnetic pressures, normalised 
here to conditions representative of the deep interiors of cool-core 
clusters. The recent observation that magnetic fluctuations in the 
solar wind are bounded by the firehose and mirror stability thresh- 
olds strongly supports the expect ation that equation ([9) will be 
satisfied in a turbulen t plasma (e.g. | Kasper. Lazarus & Garvll2002l : 
iHellinger et al. 2006; Bale et al.ll2009l) . 

There are two fundamental physical ways in which marginal 



^ In thi.s paper we use the terms 'microscales' and 'microscopic' to de- 
scribe processes whose lengthscales are just above the ion Larmor radius 
and whose growth rates are just below the ion cyclotron frequency. 
^ In a turbulent plasma, there will be regions of pressure anisotropy with 
both signs and so one might expect some average value of between 0.25 
and 1. See also .Sharma et aL 12006.) for further refinements of this mod- 
elling. 
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stability to microscale instabilities can be maintained in a weakly 
collisional plasma. Either microstabilities give rise to some ef- 
fecti ve particle scat tering mechanism that isotropises the pres- 
sure ( ISchekochihin& Cowley 2006; Sharma et al. 2006, 2007) or 
they modify on the average the structure and time evolution of 
the magnetic field (or, equivalently, of the rate of strain) so as 
to cancel the pressure anisotropy caused by the changing fields. 
While there is no complete microphysical theory, existing calcula- 
tions of particular cases ( Schekochih in et al, 2008: Califano et al.. 
20081: llst omin. Pokhotelov & Balikhirj | ^ 



200S 



Rincon. Schekochihin & Cowlev. .2010, in preparation !) all suggest 
the latter scenario. In this work we will assume that this is what 
happens. Importantly, this means that the collision rate is not mod- 
ified by the instabilities. 

Let us now use this assumption to calculate the heating rate 
(per unit volume) due to parallel viscous dissipation of motions. 
From equations l|6j and l[8j, this rate is given by 



-cTi: Vu = piAi 



bb-.Vu — ^V-u 



(11) 



We seek to write this equation solely in terms of the ion pressure 
anisotropy, and therefore we require an additional equation relat- 
ing the rate of strain (the term in parentheses) to the ion pressure 
anisotropy. This is provided bv lBraginskiil ( ll965l) : 



i^iiAi = 2.9 ( bb-.Vu - -V-u 



(12) 



This equation states that the rate of strain due to turbulent mo- 
tions generates a pressure anisotropy that is relaxed on the ion- 
ion collision timescale. We interpret the right-hand side of equa- 
tion l lI2b as including contributions to the (average) turbulent rate 
of strain from b oth the microscopic and macroscopic motions (see 
iRosin et al 1 1201 for a detailed kinetic calculation of one example 
where this is justified). Under the assumption that the collision fre- 
quency is not modified by the instabilities, equation il2l . together 
with equation states that microscale fluctuations adjust them- 
selves in response to large-scale fluid motions in order to maintain 
and preserve marginality. 

Using equation l lI2t in equation i ll It , we find that 



0+ = 0.35pi//iiAf. 



(13) 



In other words, the ion pressure anisotropy is a source of free 
energy that is eventually converted into heat by collisions. With 
marginality to microscale instabilities imposed according to equa- 
tion equation l |13l l becomes 



Q+ = 0.35 pi/^ii 



2.2 X 10"^ C^B* — . 

Pi 



(14) 



Normalised to conditions representative of the centres of cool-core 
clusters, this parallel viscous heating rate (per unit volume) is 



B 



10 fiG 



T 



2keV 



5/2 



-1 -3 

erg s cm 



(15) 



Note the lack of an explicit dependence on density and also, re- 
markably, on either the amplitude or the rate of strain of the tur- 
bulence. However, the strong dependence on B does encode infor- 
mation about all of these. As the turbulent velocity increases, so 
too will the strength of the (dynamo-generated) magnetic field (see 
Section [J!2l (, and the dissipation rate (eg. 115b will increase accord- 
ingly. In other words, the dissipation rate can be self-regulating. 
This local self-regulation is completely independent of whatever 



global self-regulation may be offered by the external source that is 
providing the turbulent energy (e.g. the AGNs), as long as there is 
enough turbulence to maintain the pressure anisotropy at its stabil- 
ity threshold (eq.|9}. 

It should be stressed that the physical picture encapsulated 
by equation il5l is somewhat different from the conventional ap- 
proach to turbulence theory. Namely, in standard turbulence theory, 
one assumes that the energy input from external driving is fixed (i.e. 
given) and that all of that power is converted into turbulence and ul- 
timately microphysically dissipated. In other words, the dissipation 
rate will adjust according to the rate at which external energy is 
input. In the ICM, there is a microscale constraint on the total tur- 
bulent rate of strain, namely, equation l|9j. Although the dissipation 
can self-regulate via the strength of the magnetic field as described 
above, not all of the power provided by the external driving in fact 
need be locally thermalised via turbulence. Indeed, the turbulence 
may have an effective 'impedance' and only accept the amount of 
power that can be locally viscously dissipated without triggering 
the microinstabilities. A further discussion of these points is given 
at the end of Section|4] 

2.3 Thermal balance and stability 

The similarity between the coefficients in equations lO and il5\ is 
striking and strongly suggests that parallel viscous heating can off- 
set radiative losses in the deep interiors of cluster cores. However, 
it is not enough simply to find a heating source capable of balanc- 
ing cooling. One must also ensure that the resulting equilibrium 
is thermally stable. This can be shown for our proposed heating 
mechanism as follows. 

In Fig.[T] we show the rates of parallel viscous heating (solid 
line) and radiative cooling (dashed line) at fixed magnetic and 
thermal pressures as a function of temperature. The thermal equi- 
librium point where heating balances cooling is denoted by a 
black dot. The cooling rate (dashed line) scales as T~''^^^ in the 
Bremsstrahlung regime, while the parallel viscous heating rate 
(solid line) scales as Hence, if the temperature is perturbed 

downwards, the heating rate increases faster than the rate of radia- 
tive cooling, and net heating restores the temperature to its equi- 
librium value. On the other hand, if the temperature is perturbed 
upwards, the heating rate decreases faster than the rate of radiative 
cooling, and net cooling restores the temperature to its equilibrium 
value. In other words, parallel viscosity, regulated by the growth 
of microscale instabilities, endows the large-scale plasma with a 
source of viscous heating that makes the plasma thermally stable. 

This is in stark contrast to a thermal balance between ra- 
diative cooling and thermal conduction (dotted line), whose heat- 
ing rate scales as T''/'^ at fixed thermal pressure. The result- 
ing equilibrium is unstable, with an isobaric perturbation result- 
ing iri_eithei_an_jsothem profile (e.g. see fig. 7 
of IParrish. Ouataert & Sharmal I2OO9I) or a cooling catastrophe. 
Thus, Spitzer conduction by itself cannot balance radiative cool- 
in g in a stable way (e.g. Bregman & David 1988; Soker 20q3; 
[Kim & Naravaij|2003bh . 

The thermal stability argument for a balance between parallel 
viscous heating and radiative cooling may be refined by allowing 
for the possibility that the magnetic-field strength is not fixed, but 
rather depends on density and temperature in some unknown way. 
(For the remainder of this Section, we suppress the subscript i on 
n and /3 for economy of notation.) Again consider a thermal equi- 
librium, denoted by a subscript '0', in which energy losses balance 
energy gains: 
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equilibrium 
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Figure 1. A sketch of the rates of parallel viscous heating (solid line) and 
radiative cooling (dashed line) versus temperature T at a fixed pressure. The 
thermal equilibrium point where heating balances cooling is denoted by a 
black dot. For a positive temperature perturbation from equihbrium, radia- 
tive cooling dominates over parallel viscous heating, whereas for a negative 
temperature perturbation from equilibrium, parallel viscous heating domi- 
nates over radiative cooling. Hence the equilibrium is stable to isobaric per- 
turbations. By contrast, thermal balance between thermal conduction (dot- 
ted line) and radiative cooling is unstable. 



Q-{no,To)^Q+{To,Bo) = Qo, (16) 

where is given by equation ([TJ and by equation J14b . Al- 
low small perturbations to this equilibrium, denoted by 5, which 
preserve a constant total pressure P: 



SP = 5p+ + i (Sp± ~5pu) =0. 



(17) 



This effectively removes sound waves from the analysis (see eq. 
[TJ and allows us to focus on n on-propagating, subsonic ('conden- 
sation') modes (e.g. Field 1965). With marginality to microscale 
instabilities imposed (eq.[9l(, equation ( I17t becomes 

SP = 5p + C-—=0, (18) 
47r 

where = 1 + 2^/3. Put simply, the pressure anisotropy modifies 
the magnetic pressure. Then the perturbed cooling and heating rates 
are, respectively. 
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If we define 
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(19) 



(20) 



(21) 



we find that the fractional perturbation in the density is related to 
the fractional perturbation in the temperature by 

5n _ ST 1 + 2C9t/Po 
no ^ To l-f 2C6'„//3o ' 

Note that this reduces to requiring a constant thermal pressure in 
the limit Po ^ 1. Substituting equation \22l into equations il9\ 
and ll20t, we find 



(22) 



SQ- _ (^T / 1 1 + 2C6It//3o 
Qo p To\2 l + 2Ce„//3o 



(23) 
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(24) 



From this, it follows that the loss function (i.e. heat losses minus 
heat gains) is 



£{n,T,B)\j 
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Since /3 > 1 in the ICM (see eq.fTOt, 

ST " 



C{n,T,B)\p ~ Qo 



To 



1 - 4 St 



(26) 



To ensure thermal stability, the term in square brackets must be 
non-negative, so that the loss function has the same sign as ST. 
Then, by equation l l26t . 



dlnB 



dlnT 



dlnB 



din 1 



< 



1 



(for thermal stability). 



(27) 



This inequality is marginally satisfied (by definition) if heating bal- 
ances cooling. Unless the perturbed magnetic field is a strongly 
increasing function of temperature or the magnetic field increases 
with decreasing density, neither of which are particularly likely, the 
inequality ([27} is satisfied and so the equilibrium is thermally sta- 
ble. 



3 IMPLICATIONS: CLUSTER EQUILIBRIUM 
PROFILES 

Given the apparent ability of parallel viscous heating to balance ra- 
diative cooling and stabilise the centres of cool-core clusters, it is 
tempting to ascribe the entire temperature profile inside the cool- 
ing radius to a long-term balance between this form of heating 
and cooling. In this Section, we investigate what this implies for 
magnetic-field strengths, turbulent velocities and turbulence scales 
throughout cluster cores as functions of radius r. We then assess a 
posteriori whether or not the results are observationally permissible 
and theoretically sensible. 



3.1 Magnetic fields 

We assume that parallel viscous heating (eq. I15t due to turbulent 
dissipation balances radiative cooling (eq. O at all radii inside the 
cluster core: 



Q+{r)~Q-{r). 

In the Bremsstrahlung regime, this implies 

\ 3/4 



B ~ 11 



^ VO.lcm-3/ V2keV/ 



(28) 



(29) 



Note that the predicted magnetic-field strength is thus a function of 
density and tem perature. It is comm on in both numerical modelling 
of clusters (e.g. Dolag et al. '200 1") and in data analysis ai ming to 
reconstruct magnetic-field str engths and spectra (e.g. Murg ia et al.l 
l2004lKuchar & Enfilin|[2()0^ to assume an exclusive relationship 
between B and no. Our arguments suggest that these models may 
need to be generalised to accommodate the temperature depen- 
dence. 
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Table 1. Predicted central magnetic-field strengths i?c, theory for ^ variety 
of clusters. Central electron number de nsities ne.c and te mperatures Tc 
for ea ch cluster are taken from (A1835) Sanders et al. (2010); (Hydra A) 
David et al. 12001); (A478, A1795) Dennis & Chandran (2005); (A2199) 
Johnstone et al. (2002); (M87) lOhizzardi et al. (200#; (Cent a urus, A262) 
Cavagnolo et aU ( .2009,) ; (A2142, A 401) .MarkevitcM h998h. IChen et al, 
j2007t) ; (Opiiiuchu s, A2634, A4 Q) IChen et al.l <2007ir Fukazawa et al 
119981) ; (A2382) lEbehng et alj il996l) iGuidettietall feOOSi); ari d 
(A2255) lFeretti et alj h997h . lGovoni eiM j2006l) . iDavis & White! ( Il998l) . 



Observationally-inferred magnetic-field strengths Sc.obs. vi'here available, 
are taken from (Hydra A) [ Vogt & EnBlin (20()1; (A2199. M87, A1795 



A2634 , A400) Eilek & Owen (2002); (Centaurus) "Xavlor. Fabian & Allei 
i2002l) ; (A2382) Guidetti et al. (2008); and (A2255) Govoni et al. (200( 
The non-cool-core cluster results are discussed in Section[33] 



lie3 



Cluster name 


"o.c 

(10-2 cm-3) 


(keV) (5- 


, theory 

1/VG) 


-^c,obs 

(mG) 


Cool-core clusters 


A1835 


10 


2.85 


13.8 




Hydra A 


7.2 


3.11 


12.4 


12" 


A478 


15.2 


1.72 


12.1 




A2199 


10 


~ 2 


~ 11 


15" 


M87 


10.8 


1.62 


9.8 


35" 


A 1795 


5.4 


2.26 


8.6 


9.7" 


Centaurus 


9.5 


1.24 


7.7 


8 


A262 


3.7 


1.54 


5.5 




Non-cool-core clusters 


A2142 


1.87 


8.8 


13.0 


RM<^ 


Ophiucus 


0.80 


10.3 


9.5 


RM<= 


A401 


0.70 


8.3 


7.6 


RM<= 


A2382 


0.50 


2.9 


3.1 


3 


A2634 


0.28 


3.7 


2.7 


3.5" 


A2255 


0.2 


3.5 


2.2 


2.5 


A400 


0.24 


2.3 


1.8 


2.9" 


Estimates for Bc.obs in Hydra 


A are widely 


varying. 


from 7 fiG 



to ~ 30 AtG (Taylor & PerlevHl993l;lKuchar & Enfihn.20091) . 
^ For these clusters, oba is inferred by assuming that only one mag- 
netic filament along the line of sight accounts for the observed rotation 
measure; therefore these are likely upper limits. 

'RM' indicates that a rotation mea sure exists, but that no Bt.,obs iias 
been inferred (see lGovoni et alj201(l and references therein). 



For typical electron densities (ric ~ 0.01 - 0.1 cm- ) and 
temperatures (T ~ 1-3 keV) at the centres of cool-core clus- 
ters, equation i29\ implies central magnetic fields ~ 1 - 10 fiG, 
within observational constraints. For example, conditions near the 
centre of the po pular Hydra A clu ster (nc,c — 0.072 cm-^ and 
Tc ~ 3.11 keV: iDavid et alj|200ll : H. Russell, private communi- 
cation) imply a thermal-equilibrium magnetic-field strength Be — 
12.4 ^^^/2 paj-tfier out around ~ 30 kpc, the observed elec- 
tron density — 0.02 cm~^ and temperature T ~ 3.5 keV 
imply i? ~ 7 /iG. These are both in good agreement 

with magnetic- field strength estimat e s in H ydra A from Faraday 
rotation maps dVogt&EnBiii] 12003 '. 2005). For anothe r cluster, 
A2199, popular among theorists (e.g. Parrish, Ouataert & Sharm^ 
|2009|) . a cent ral density rte.c — 0.1 cm-'' and central temperature 
Tc ~ 2 keV dJohnstone et"ai]|2002l : H. Russell, pr i vate c ommuni- 
cation) imply Be ^ 11 C"^^^ uG. lEilek & Owed ( |2002|) inferred 
a central magnetic-field strength there of 15 fiG by assuming that 



only one magnetic filament along the line of sight accounts for the 
observed rotation measure. In Table [T] we list these and other cen- 
tral magnetic-field strength predictions. 



3.2 Turbulent velocities 

In order to estimate the turbulent velocities in the ICM, we assume 
that the large-scale kinetic and magnetic energies are in overall 
equipartition, 

(30) 



1 rr2 



Stt' 



where Uims = (u^)^^'^ is the rms flow velocity. This is ex- 
pected to be the case for a magnetic field amplified and brought 
to saturation by the fluctuati on dynamo (see, e.g., review by 
ISchekochihin & CowlevI |2007| , and references therein). Then the 
rms turbulent velocity is equal to the Alfven speed and, using equa- 
tion i29\, we therefore obtain 



C/rr 



70 f 



-1/2 



T 



2keV 



\ 3/4 

J km s~ 



(31) 



in the Bremsstrahlung regime. In other words, relatively hotter 
(cooler) cores should therefore have larger (smaller) turbulent ve- 
locities. The corresponding Mach number is 



M ■-. 



0.18 C 



-1/2 



2keV 



1/4 



(32) 



where = (T /ms)^^"^ is the sound speed. Note the weak depen- 
dence of M on temperature. For T ~ 1-10 keV, f/rms ranges 
from ~ 44 C^^'^ km s"^ to 210 C^^'^ kms-\ so M ranges 
from ~ 0.16 ^"^^^ to 0.24 While the turbulent velocities 

are no t yet measured d i rectly , t heoretical ( e.g. Dennis & Chandrati 
2005'; lEnfilin & Vog3 
2006), numerical (e.g. 



20061: 



Subramanian, Shukurov & Hauger 
Norman & Bryan .1999.: .Ricker & Sarazin 



20()]J; ISunvaev. Norman & Brvanj l20()3h and indirect obser- 
vational (e.g. Schuecker et al. 2004; Churazov et al. 200'" 
Rebusco et al. 2005, 2006; Graham et al. 2006; Rebusco et 



20081 ; ISanders^ Fabian & Smithli2010l) estimates suggest that the 



numbers we are predicting are reasonable. 

We caution here that equation ( I3U should be considered 
a lower limit on the actual turbulent velocity since C/rms is 
unlikely to be smaller, but could be larger, than the Alfven 
speed. In magnetohydrodynamic numer ical simulations, it is of- 
ten larger by a factor of order unity (e .g. ISchekochihin et al]|2004l ; 
iHaugen, Brandenburg & Dobleill2004h . 



3.3 TXirbulence scales 

Implicit in the above discussion is the requirement that there be 
enough turbulent energy for the viscous heating rate mandated by 
the marginal stability condition (see eq. I14t to be maintained. We 
assume that throughout the cluster core all of the power from exter- 
nal stirring (see discussion in Section|4ll accepted by the turbulence 
is locally dissipated and thermalised via parallel viscosity; 



mini 



2 

rms 



Tturb 



(33) 



'turb '^he effective rate at which energy is converted 
into turbulent motions. In conventional turbulence, this is equiv- 
alent to the eddy turnover timescale. In writing equation ( I33t , we 
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Figure 2. (Le/f to right, top to bottom) Profiles of the observed (deprojected) electron number density ric, observed (deprojected) temperature T, predicted 
magnetic-field strength B, predicted rms turbulent velocity Urms, predicted turbulence lengthscale L and predicted turbulent diffusion coefficient Kturb for 
the cool-core cluster A1835. The predicted data points correspond to \^\ = 0.75, halfway between the firehose instability threshold (\^\ = 1; dotted line) 
and the mirror instability thi'eshold = 0.5; dashed line). The line plots are derived from best-fitting analytic profiles to the electron numb er density and 
tempe rature (see footnote[3). The thick arrow on the plot of Uims denotes the 182 km observationally-derived upper limit obtained bv lSanders et alj 
i2O10l) at r ~ 30 kpc. 



have ignored the possibihty that energy could cascade to conision- 
less scales via Alfvenic turbulence and heat the plasma via mi- 
crophysical dissipation a t the ion and electron Larmor scales (see 
ISchekochihin et alj|2009l and references therein). 
Equations iSll and i33\ give 

---2r^(^)-^(^)Myr. (34) 

We also define the characteristic turbulence lengthscale ('outer 
scale') 

L = Urnis "^turb 

Another corollary of equations ( 1311 ) and i34\ is that the turbulent 
diffusion coefficient is 

r-r2 

'^turb tyj-jng'^turb 

~ o in27 c-2 ( We ( T \^'^ 2 -1 

~ 3 X 10 t — ; — -- cm s 

^ VO.lcm-3/ V2koVy 

(36) 

Note that Atturb has the same scahng as the lSpitzeJfTgei) electron 
thermal diffusion coefficient (see Section [T6t . Equations jiW and 
([35} further imply an effective Reynolds number associated with 
the parallel viscosity 

(Kvisc is the viscosity coefficient; see lBraginskiJl96^ . Thus, Re ~ 
1-10 and is independent of radius. In other words, the outer and 
viscous scales are close to one another, so that the motions are 



dissipated near th e outer scale and there is no inertial range (c.f. 
Ifabian et al. 2003). However, we caution that an effective Reynolds 
number of order unity does not imply laminar flow in this case, 
since the turbulence is randomly stirred. We also stress that this is 
Re calculated on the basis of parallel collisional viscosity - it does 
not imply a viscous cutoff for all plasma fluctuations, although, as 
stated above, we have ignored those as a possible heating channel. 

3.4 Cool-core cluster profiles: A1835 

Using equations ( |29l ) - ( 1361 ), profiles of B, f/rms, L and Kturb may 
be calculated for any given cluster. Here we provide results for one 
particular cool-core cluster, A 1835 (we have performed the same 
exercise for many other clusters, with similarly sensible outcomes). 
In Fig.m we give the observed (deprojected) electron number den- 
sity Tie and the observed (deprojected) temperature T, which are 
then used to predict the magnetic-field strength B, the rms turbu- 
lent velocity (7rms, the characteristic turbulence scale L and the 
turbulent diffusi on coefficient Kturb. The density and temperature 
profiles are from lSanders et all ( I2OIOI) . The data points for the pre- 
dicted quantities use |^| — 0.75, but we also show upper and lower 
limits calculated by setting |^| = 0.5 (dashed lines) and j^j = 1 
(dotted lines), respectively. 

Equations i29l - i36l do not imply a specific causal relation- 
ship between the five quantities Uc, T, B, [/rms and L. Given any 
two quantities, our theory can predict the other three. For example, 
if observations determined the profiles of Uc and ?7rms for a given 
cluster, rather than of ric and T, then our theory would predict T, B 
and L. Ideally, our theory would be put to the test if three or more 
of these profiles were known observationally. 

The predicted value for B near the centre of the core is ~ 15 - 
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21 /iG, decreasing to ~ 5-7 at the outer core boundary. To our 
knowledge, observational estimates of the magnetic-field strength 
in A 1835 have not yet appeared in the literatu re. As a radio mini - 
halo has recently been detected in A1835 ( Murgia etalj|2009l) . 
there is hope for a magnetic field measurement there, although 
this might be quite a difficult task for such a distant cluster. The 
predicted turbulent velocity dispersion Kms ~ 70 - 270 km 
throughout the core, attaining a value C/rms — 114 - 162 km at 
a radius of 30 kpc. Thi s is wi thin the 182 km upper limit ob- 
tained bv lSanders et alj j201Cl) by measuring emission lines within 
r ~ 30 kpc, which is denoted on the plot by a thick arrow. A 
more conservative observational estimate of f/rms ^ 274 km s^^ 
was derived by treating the cl uster as a poin t source and not ap- 
plying any spatial smoothing jSanders et aL|[20im . The predicted 
turbulence scale isL ~ 0.2-0.7 kpc near the centre of the 
core, increasing outwards to ~ 70 - 200 kpc near the temperature 
maximum. The predicted diffusion coefficient Kturb rises sharply 
from ~ 10 cm s near the core centre to ~ 10 cm s at 
the outer core boundary. Given the uncertainties discussed in Sec- 
tions [3]2] and [53] it is encouraging that these values are compara- 



in a variety of observe d clusters (Rebusco et al.ll200^l2006ll2008l : 
iDavid & Nul"sedll2008l) . 



3.5 Non-cool-core cluster profiles 

The plasma microphysics responsible for parallel viscous heating 
is in no way unique to cool-core clusters. Since plenty of turbu- 
lence is expected to be present (stirred by mergers, etc.), pressure 
anisotropics will develop and will presumably be maintained at a 
marginal level as described in Section lZ2l While heating and cool- 
ing times are quite long in such clusters and the situation can be 
quite time-dependent, we would nevertheless like to explore what 
the conjecture of an approximate local heating-cooling balance 
would imply, with the caveat that such a balance may in princi- 
ple take a very long time to be established. We will see that such a 
balance leads to quite reasonable predictions that seem to be borne 
out by observational data. 

In isothermal clusters, equations {19\ - ( |35t imply that 

D 1/2 



const, 



Tturb OC L cx rip 



(38) 
(39) 
(40) 



The scaling B oc n^^'^ in non-cool-core cluste rs has in fact 
alread y been observationally infer red for A2 382 ijGuidetti et al.l 
|2008|) . Coma jSonafede etallbOlOl) and A665 l lVacca et alj|2010l) . 

The central density n e.c — 5 x 10~^ cm ~^ and temperature 
T ~ 2.9 keV in A2382 jEbeling et ai]| 19961) implies a thermal- 
equilibrium magnetic-field strength Be — 3.1 ^"^/^ /iG, in excel- 
lent agreement with the_2^_3jiG estimate derived from rotation 
measure observations jOuidett i et al. 200^ For A2255, rie.c — 



2 X 10"^ cm^fjFeretti et al. 19971 : iGovoni et al. 2006) and T 



3.5 k eV jPavis & White. 19981 however, see Sakelliou & Ponmanl 
I2OO6I who find temperature variations across A2255 from T ~ 5.5 
- 8.5 keV), the thermal-equilibrium magnetic-field strength Be — 
2.2 ^"^/^ ^G. This compares fav ourabl y with the observational 
estimate B ~ 2.5 /iG obtained bv lGovoni et al.l (|2006j), who used 
B oc in their analysis. In Table [T] we list these and other 

central magnetic-field strength predictions. 

Very recently there have been observational estimates of the 



m agnetic-field strength p rofile in the non-cool-core Coma cluster 
bv lSonafede et alj j2010l) . They found that the best-fitting profile 
of the magnetic-field strength is 



4.7 



3.44 X 10-3 cm-3 



1/2 



mG, 



(41) 



with a '/3-moder jCavaliere & Fusco-Femianoll 19761) electron den- 
sity profile 



no = no 1 + ^ 



-3/3/2 



(42) 



where no = 3.44 x 10"^ cm"^, Tc = 291 kpc and /3 = 0.75. 
Taking Coma to be an isothermal cluster with temperature 8.2 keV 
dAmaud et al.l200]l) . we find that the implied parallel viscous heat- 
ing rate is (from eq.[T4]using eq. l41b 



)+ ^ 1.6 X 10' 



3.44 X 10-3 cm- 



erg s cm . (43) 



By way of comparison, the implied Bremsstrahlung cooling rate is 
(from eq.[3j 



Q- = 2.5 X 10"^* 



ne 



^3.44 X 10-3 cm- 



erg s ^ cm , (44) 



where Uc is given by equation i42t . While it is rather curious that 
these two rates are not only comparable, but also have the same 
radial scaling, we are not on solid ground applying our estimates to 
Coma because the cooling time for it implied by equation l l44t is 
icooi ~ nT / Q- ~ the Hubble time at the cluster centre, increasing 
outwards in radius. The other isothermal clusters mentioned above 
and in Table[T]have relatively shorter cooling times, so one could in 
principle have a heating-cooling balance - but no radial dependence 
of the magnetic-field strength has so far been measured. 



3.6 Postscript: thermal conduction 

Bearing in mind that thermal conduction models routinely 
fail in the innermost region s of coo l cluster core s (e.g. 
Zakamska & Narayan 2003; Mar kevitch et al. 2003; Kaast ra et al.l 
I2OO4I ; iGhizzardiet all 12004 ). it is important to note that parallel 
viscous heating should be especially important in these relatively 
cold (T ~ 1 keV) and strongly-magnetised [B ~ 10 /iG) regions. 
A balance between parallel viscous heating and radiative cooling, 
however, does contain the implicit assumption that the thermal con- 
duction i s relatively un important. This ought to be checked. 

The lSpitzeJ ( Il962i) electron thermal diffusion coefficient is 



«:Sp 



10^ 



0.1 cm- 



r 



2keV 



5/2 



(45) 



In tangled turbulent magnetic fields, the true thermal conductivity is 
expected to be k = /ksp^ where / is the suppression facto r, which 
is expected to ran g e anywhere between ~ 10~^ and ~ 0.3 (Tribbld 



'Kol I1995I: [Pistinner & Shayivl Il996l: IChandran & Cowle 



1998; Malvshkin & Kulsruj '2OOII; iNaravan & MedvedevI l200lF 
Gruzinov 2002; Cho et al. 2003). Note that the Spitzer diffusion 
coefficient has the same density and temperature scaling as the tur- 
bulent diffusion coefficient derived in Section [33] We see that for 
/ < 0.03 turbulent heat diffusion is likely to be at least as 
important as the coUisional conductivity. 
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Figure 3. (a) The cooling/heating timescale tcool/hcat (solid line), ion-electron equilibration timescale ti_c,cq multiplied by 10"^ (dash-dot line) and the 
shortest possible conduction timescale tcond.min (dashed line) for A 1835. (b) The Spitzer conduction suppression factor / = k/ksp required for thermal 
conduction to balance radiative cooling as a function of radius for several clusters. Only suppression factors / < 1/3 (shaded region) are allowed in the 
presence of a tangled magnetic field. The line plots are derived from best-fitting analytic profiles to the electron number densities and temperatures (see 
footnote [3). 



In Fig. [3^, we plot for A 1835 the cooling/heating timescale 



^cool/hcat 



3:^ = 0.21 
2Q- Vo.li 



T 



2keV 



1/2 



the ion-electron temperature equilibration timescale 



9.3 



( V 

VO.l cm-»/ 



T 



2keV 



3/2 



kyr, 



and the shortest possible conduction timescale 



^coiidjinin — 



}_d_ 

^2 Q^. 



fUcKSpr 



dr 



Gyr, 
(46) 

(47) 
(48) 



that one can expect in the presence of a t angled magnetic field 
(with / = 1/3; see, e.g.. iGruzino Heating via thermal 

conduction is everywhere unimportant relative to radiative cool- 
ing and (by construction) parallel viscous heating. Also note that 
ti_e,eq ^ 1 Myr everywhere in the cluster core, so that the as- 
sumption of equal ion and electron temperatures is well justified. 

These conclusions are by no means unique to A 1835. To 
demonstrate this point, we plot in Fig. [3}3 the Spitzer conduction 
suppression factor / required for thermal conduction to balance ra- 
diative cooling as a function of radius for several different clusters[f| 
Only suppression factors / < 1/3 are expected in the presence of 



^ For the line plots in the figures, we have used analytic fits to the observed 
electron density and temperature profiles. One advantage is that this en- 
sures smooth gradients for computing the c onductive heating rates. The pro- 
files for A478 and A1795 were taken from Dennis & Chandran I 200|) and 
the profiles for M87 were taken f rom iGhizzardi et a l. ( 2004). For A1835 
iSanders et ai]|2010t) and Hydra A jPavid et alBoOlh , the electron number 
density was fit using equation )42t with /3 = 0.593, = 32.42 kpc and 
no = 0.115 cm-3 (for A1835) and /3 = 0.393, Vc = 10.9 kpc and 
no = 0.0669 cm~^ (for Hydra A). Their temp erature profiles were fit 
using equation (22) of IPennis & Chandraij j2005h with To = 9.55 keV, 
Ti = 7.39 keV, ret = 33.35 kpc and 5 = 0.62 (for A1835) and 
To = 3.95 keV, Ti = 0.856 keV, m = 30.60 kpc and (5 = 0.48 
(for Hydra A). 



a tangled magnetic field. None of these clusters can be in thermal 
balance between thermal conduction and radiative coolingQ 



4 DISCUSSION 

In this paper, we have introduced a model for regulating cooling 
in cluster cores in which turbulence, magnetic fields and plasma 
physics all play crucial roles. Our findings are fundamentally based 
on an appreciation that, whatever the source of effective viscos- 
ity in the ICM, it is certainly not hydrodynamic. Instead, it is set 
by the microscale plasma-physical processes, which are inevitable 
in a weakly coUisional, magnetised environment such as the ICM. 
Assuming that microscale plasma instabilities pin the pressure 
anisotropy at its marginally stable value, we derive an expression 
for the heating rate due to parallel viscous dissipation of turbulent 
motions. For typical conditions in a variety of cluster cores, this rate 
is of the right magnitude to balance radiative cooling. Moreover, 
this source of heating turns out to be thermally stable. Put simply, 
the viscous stress in the ICM is a dynamic quantity that responds to 
local changes in temperature, density and magnetic-field strength in 
such a way as to prevent runaway heating or cooling. A basic qual- 
itative outline of how this occurs is given in Fig. |4l If true, what 
we have conjectured constitutes a physical mechanism that allows 
clusters to develop stable non-isothermal temperature profiles and 
avoid a cooling catastrophe - an outcome that has remained elusive 
in models involving balancing the cooling by thermal conduction. 

Of the five quantities - density, temperature, rms magnetic- 
field strength, rms turbulent velocity, characteristic turbulence scale 
- given radial profiles of any two, we can predict the rest. The 
most reliable, as well as readily observable, predictions that follow 
from our theory concern the (rms) magnetic-field strength. For typ- 
ical electron densities and temperatures, we predict magnetic-field 



* Note that there do exist some clusters for which a fraction of Spitzer con- 
duct ivity may be a ble to stabihse their cores against a coohng catastro phe 



o pne 

(e.g. Zakamska & Naravan 2003; see however § 3.3 of .Conrov & Ostriked 
I2OO8 for an alternative point of view.) 
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energy source (AGNs, mergers, etc.) 



energy absorption rate (effective impedance) 



C^ Aulence: U^^ ^ 




pressure anisotropy pinned at marginal stability: ~ /^^ ' 




thermal balance 



Figure 4. Qualitative outline of the processes responsible for thermal bal- 
ance between parallel viscous heating and radiative cooling. Energy sources 
(AGNs, mergers, etc.) inject energy into the ICM, some or all of which (de- 
pending on the effective impedance of the ICM) is absorbed by the plasma 
and converted into turbulence. This changes the magnetic-field strength, 
giving rise to both a fluctuation dynamo and pressure anisotropics. The fluc- 
tuation dynamo presumably saturates in equipartition between the turbulent 
magnetic and kinetic energies. The pressure anisotropics excite microscale 
instabilities, whose effect is to pin the pressure anisotropy at marginal sta- 
bility. The pressure anisotropy determines the viscous stress and therefore 
the heating rate. This heating balances radiative cooling, giving rise to a 
stable thermal equilibrium that is maintained by energy injection by self- 
regulated sources (such as AGNs) and/or self-regulated energy absorption 
by the turbulence. 



strengths in the range --^1-10 /iG, well within current obser- 
vational constraints. For the specific clusters discussed here (e.g. 
A1835, Hydra A, A2199, A2382, A2255), a balance between par- 
allel viscous heating and radiative cooling results in field strengths 
and profiles that are quite reasonable and in good agreement with 
current observational estimates where available. It would be inter- 
esting to test our model further via analysis of Faraday rotation 
maps. 

Another prediction is the manner in which the magnetic-field 
strength B depends on electron density rio and temperature T (eq. 



B oc nl'^T'"^ 



(for cool-core clusters). 



(49) 



In cool-core clusters, this suggests that the oft-employed assump- 
tion of an exclusive relationship between B and rie should be re- 
laxed. It is encouraging that, in isothermal clusters, the implied 
scaling 

D 1/2 

B oc n' 



(for isothermal clusters) 



(50) 



has already been o bservation ally inferred to hold in A2382 
tauidetti etal.' '2008'), Coma teonafede etall |2010|) and A665 
(Vaccaetal. 2010). A caveat here is that a local heating-cooling 
balance in such clusters may take very long to establish and so 
should be treated with due scepticism. A test of whether this con- 
jecture is reasonable would be to determine whether equation ( |50t 
is satisfied systematically in a large sample of isothermal clusters. 



There are several peripheral consequences of our predic- 
tions that warrant further discussion. Firstly, the relatively large 
magnetic-field strengths inferred from observed rotation measures 
and predicted by a balance between parallel viscous heating and 
radiative cooling are stable or at le ast marginally s table to the HBI. 
The stability criterion for the HBI ( IOuataerlll2008l) may be written 
as a lower-bound on the magnetic-field strength: 



S > 3 



1/2 



10"*cms-2/ VO.lcm-3 
AlnT/AlnrV''^ ^ 

" ' 0.1 



1/2 



100 kpc 



1/2 



(51) 



Even under the favourable conditions we have chosen for the elec- 
tron density nc, the radial distance from the cluster centre r and the 
temperature slope dlnT / d\nr, the magnetic fields predicted in 
this pa per exceed the stability limit (15 lb. Indee d. Bo gdanovic et al.l 
l l2009l § 3.3) and lParrish, Ouataert & Sharmal ( l2009l g 5.7) found 
that simulated clusters with magnetic-field strengths B ~ 3 
demonstrated a delayed cooling catastrophe due to the stabilisation 
of HBI modes by magne tic tension. By contrast, the rms ma gnetic- 
field strength used in the lParrish. Ouataert & Sharmal ( 1201 Oh simu- 
lation, where a combination of turbulent stirring and HBI-govemed 
thermal conduction gave a cool core in long-term thermal bal- 
ance, was 10"" G. More observational estimates of magnetic-field 
strengths in a variety of clusters would clearly be very useful to 
help realistic modelling and testing of theories. 

Secondly, the same process that is responsible for stably heat- 
ing the ICM in our model may also influence the turbulent diffu- 
sion of metals throughout the cores of galaxy clusters. That our 
values for the turbulent diffusion coefficient Ktmb are comparable 
to those inferred in a variety of clusters is encouraging. However, 
iRebusco et alj ( |2005|) found that diffusion coefficients that increase 
with radius (as ours does) imply abundance profiles that are too 
centrally peaked compared to those observed. It is therefore worth 
investigating the diffusion of metals using our predicted scaling 
Kturb OC n-^T^/^ 

Thirdly, turbulent heat diffusion may become dominant rela- 
tive to collisional electron heat conduction if the Spitzer conductiv- 
ity suppression factor / < 0.03 (see eansl36landl45t. However, for 
typical density and temperature profiles of cool-core clusters, nei- 
ther form of heat diffusion seems to be important relative to parallel 
viscous heating (see Section lT6t . 

It is prudent to reiterate here our assumptions and their limita- 
tions: 

(i) We have implicitly assumed that there is enough turbulent 
energy so that, if it is thermalised via parallel viscous heating, it 
can offset cooling. This requires either the source of the turbulent 
energy or the amount of energy locally accepted by the plasma and 
converted into turbulent motions to have some knowledge of the 
cooling rate and to be self-regulating. 

AGNs are a natural candidate for providing a self-regulating 
energy source, whether the stirring is due to AGN-driven jets, 
bubbles, weak shocks/sound waves, gravitational modes and/or 
cosmic-ray-buoyancy instabilities. Observationally, a large major- 
ity of cool-core clusters harbour radio sources at their centres, and 
the AGN energy output inferred from radio-emitting plasma out- 
flows and caviti es is often similar t o the X-ray cooling rate of the 
central gas (e.g. lBlantonetalj20I0l) . 

However, it may not be necessary for the energy source to be 
fine-tuned or tightly self-regulating. Our scheme for heating the 
ICM can deal with excess turbulent energy via the following two 
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considerations. First, not all of the power provided by the exter- 
nal stirring has to be locally thermalised via turbulence, as the tur- 
bulence may have an effective 'impedance' and only accept the 
amount of power that can be locally viscously dissipated with- 
out triggering the microinstabilities; the remaining power could be 
transported elsewhere. Second, the parallel viscous heating rate (eq. 
113b . coupled with the assumption of equipartition kinetic and mag- 
netic energies (eq. I30t , is naturally self-regulating: any increase in 
turbulent energy implies an increase in magnetic energy, which im- 
plies an increased viscous dissipation rate. Both of these regulation 
mechanisms require there to be a sufficient amount of turbulence 
to pin the pressure anisotropy at its marginal stability threshold. 
As it happens, t here seems to be no dearth of turbulent power (e.g. 
IChurazov et al.i2004) . The fact that some of this power can be ther- 
malised via parallel viscous dissipation in a thermally stable way 
provides an attractive alternative to other heating models that suffer 
from stability issues. Further work on the excitation and replenish- 
ment of turbulence in a weakly collisional ICM is clearly needed. 

(ii) The heating is assumed to be all due to parallel collisional 
viscosity. In principle, energy could cascade through the parallel 
viscous scale and onwards to collisionless scales via Alfvenic tur- 
bulence and then heat the plasma via m icrophysical dissipation at 
the ion and electron Larmor scales (see lSchekochihin et alj|200^ . 
and references therein). This possibility has been ignored here. 

(iii) The predicted turbulent velocity, obtained by assuming ap- 
proximate equipartition between kinetic and magnetic energies, is 
probably a lower limit on the actual turbulent velocity since f/rms is 
unlikely to be smaller than the Alfven speed. Unfortunately, a more 
exact estimate of t/rms would require a detailed understanding of 
turbulent dynamo saturation in the ICM, far beyond the scope of 
the present work. 

(iv) Our estimates of the characteristic turbulence scale and tur- 
bulent diffusion coefficient are perhaps the more uncertain of our 
predictions, as they depend on the rate of transfer of energy from 
the external driving sources into the turbulence. This encodes what 
we have referred to above as the effective 'impedance' of the tur- 
bulence, the detailed physics of which is poorly understood. 

All these concerns clearly require further theoretical work. How- 
ever, the fact that our predictions are not too far from current ob- 
servational estimates or plausible expectations lends us hope that 
these concerns might be of secondary importance. 

While the exact numbers predicted in this paper should be 
taken with a grain of salt, one cannot help being encouraged by the 
fact that they seem to be quite reasonable without any fine-tuning 
of adjustable prefactors. The basic conjecture that the saturation of 
microscale plasma instabilities endows the ICM with a thermally 
stable source of viscous heating seems robust. If our predictions are 
confirmed, this would constitute strong evidence that microphysi- 
cal plasma processes play a decisive role in setting the large-scale 
structure and evolution of galaxy clusters. 
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